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Abstract 

We discuss in detail solar neutrino oscillations with 5m 2 /E in the range 
[IO -10 , IO -7 ] eV 2 /MeV. In this range, which interpolates smoothly be- 
tween the so-called "just-so" and "Mikheyev-Smirnov-Wolfenstein" oscillation 
regimes, neutrino flavor transitions are increasingly affected by matter effects 
as 5m 2 /E increases. As a consequence, the usuai vacuum approximation has 
to be improved through the matter-induced corrections, leading to a "quasi- 
vacuum" oscillation regime. We perform accurate numerical calculations of 
such corrections, using both the true solar density profile and its exponen- 
tial approximation. Matter effects are shown to be somewhat overestimated 
in the latter case. We also discuss the role of Earth crossing and of energy 
smearing. Prescriptions are given to implement the leading corrections in the 
quasi-vacuum oscillation range. Finally, the results are applied to a global 
analysis of solar v data in a three-flavor framework. 

PACS number(s): 26.65.+t, 14.60.Pq 
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I. INTRODUCTION 



A well-known explanation of the solar v e flux deficit []TJ is provided by flavor oscillations 
of neutrinos along their way from the Sun (©) to the Earth (©). For two active neutrino 
states [say, (i/ e , z/ M ) in the flavor basis and (i/j, z/ 2 ) in the mass basis], the physics of solar v 
oscillations is governed, at any given energy E, by the mass-mixing parameters 5m 2 and lo 
in vacuum,[] as well as by the electron density profile N e (x) in matter ||. 

Different oscillation regimes can be identified in terms of three characteristics lengths, 
namely, the astronomical unit 

L = 1.496 x IO 8 km , (1) 

the oscillation length in vacuum 



AttE , / 5m 2 /E \ 



-1 



— — = 2.48 x 10 — 5-i km , (2) 

óm 2 \eV 2 /MeVJ V ; 

and the refraction length in matter 



-i 



£ ma t = = 1-62 x IO 4 ( — ) km , (3) 

V2G F N e ^mol/cm 3 / v ; 

which is associated to the effective mixing angle uj m 0, 

• 2 n sin 2 2cu 

sin 2uj m = — 2 ■ (4) 

(cos2u; - L osc /L mat ) 2 + sin 2lj 

Typical solutions to the solar neutrino problem (see, e.g., [[|) involve values of L osc either 
in the so-called "just-so" (JS) oscillation regime ||, characterized by 

^osc ~ L > L mat , (5) 
or in the "Mykheyev-Smirnov-Wolfenstein" (MSW) oscillation regime 0, characterized by 

L™ W - L mat « L . (6) 

The two regimes correspond roughly to 8m 2 / E ~ O(10 -11 ) eV 2 /MeV and to 8m 2 / E > IO -7 
eV 2 /MeV, respectively. 

For just-so oscillations, since L mat /L osc — > 0, the effect of matter is basically to suppress 
the oscillation amplitude both in the Sun and in the Earth (sin 2 2u m — >• 0), so that (coherent) 
flavor oscillations take place only in vacuum, starting from the Sun surface . Conversely, 
for MSW oscillations, L osc ~ L mat and flavor transitions are dominated by the detailed 
matter density profile, while the many oscillation cycles in vacuum (L osc L) are responsible 
for complete v decoherence at the Earth, once smearing effects are taken into account M. 



1 One can take dm 2 = mi, — m 2 > 0, as far as lo is taken in the first quadrant [0, |]. 
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Therefore, it is intuitively clear that in the intermediate range between ([5]) and (^j), 
corresponding approximately to 10~ 10 ^$ 5m? / E < IO -7 eV 2 /MeV, the simple vacuum 
oscillation picture of the JS regime becomes increasingly decoherent and affected by matter 
effects for increasing values of Sm 2 /E, leading to a hybrid regime that might be called of 
"quasi-vacuum" (QV) oscillations, characterized by 



v mat ~ 



< L QV < L (7) 

~ -^osc ~ J -' ) V 1 ) 



In the past, quasi-vacuum effects on the oscillation amplitude and phase have been ex- 
plicitly considered only in relatively few papers (see, e.g., as compared with the vast 

literature on solar neutrino oscillations, essentially because typical fits to solar v rates al- 
lowed only marginai solutions in the range where QV effects are relevant. However, more 
recent analyses appear to extend the former ranges of the JS solutions upwards |15| and of 



the MSW solutions downwards |TTJ in 5m /E, making them eventually merge in the QV 
range [|17]], especially under generous assumptions about the experimental or theoretical v 
fìux uncertainties.0 Therefore, a fresh look at QV corrections seems warranted. Recently, 
a semianalytical approximation improving the familiar just-so formula in the QV regime 



was discussed in W5§\ and, in more detail, in 111911, where additional numerical checks were 



performed. In this work we revisit the whole topic, by performing accurate numerical cal- 
culations which include the exact density profile in the Sun and in the Earth, within the 
reference mass-mixing ranges 5m 2 / E e [10~ 10 , IO" 7 ] eV 2 /MeV and tan 2 u; e [10~ 3 , 10] .0 We 
also discuss some approximations that can simplify the computing task in present applica- 
tions. We then apply such calculations to a global analysis of solar neutrino data in the 
range 5m 2 < IO -8 eV 2 . 

Our paper is structured as follows. The basic notation and the numerical techniques 
used in the calculations are introduced in Sec. |I| and |T|, respectively. The effects of solar 
matter in the quasi-vacuum oscillation regime are discussed in Sec. [TV], where the results for 
true and exponential density profiles are compared. Earth matter effects are described in 
Sec. [V]. The decoherence of oscillations induced by energy (and time) integration is discussed 
in Sec. |VT| The basic results are summarized and organized in Sec. |VH , and then applied 



to a three-flavor oscillation analysis in Sec. Vili. Section IX concludes our work 



II. NOTATION 

The v propagation from the Sun core to the detector at the Earth can be interpreted as 
a "doublé slit experiment," where the originai v e can take two paths, corresponding to the 
intermediate transitions v e —>■ v\ and to v e — > v 2 - The global v e survival amplitude is then 
the sum of the amplitudes along the two paths, 



2 Such extended range for current solutions reflects, in part, the lack of a strong, model-independent 
signature of energy dependence in the solar neutrino deficit. 

The variable tan uj is useful to chart the first two octants of the mixing angle range 
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A(v e -> v e ) = A Q (v e -> z/i) • A vac (z/i -> i/i) • A ffi (i/i -> i/ e ) 

+ A Q (v e -> z/ 2 ) • A vac (z/ 2 -> i/ 2 ) • A e (z/ 2 -> i/ e ) , (8) 

where the transition amplitudes from the Sun production point to its surface (A@), from 
the Sun surface to the Earth surface (A vajC ) an d from the Earth surface to the detector (A & ) 
have been explicitly factorized. The v e survival probability P ee is then given by 

P ee =\A(v e ^v e )\ 2 . (9) 

In general, the above amplitudes can be written as 

A & (v e -> v x ) = y^exp^©) , (IOa) 

A yac {u 1 ^u 1 )=exp{-iml(L-R G )/2É) , (lOb) 

A & (yi -> z/ e ) = y^exp(^ e ) , (lOc) 
for the first path and as 



A Q (v e ^v 2 ) = ^l-P Q , (Ila) 
A vac {v 2 ^v 2 )=exp(-im 2 2 {L-R G )/2E) , (llb) 



A e (^ 2 -n> e ) = V!- p e , (He) 

for the second path, where R & is the Sun radiusfj In the above equations, P Q and P ffi are 
real numbers (e [0,1]) representing the transition probability P(v e <-> Vi) along the two 
partial paths inside the Sun (up to its surface) and inside the Earth (up to the detector). 
The corresponding phase differences between the two paths, £ and £ e (e [0, 2%]), have been 
associated to the first path without loss of generality. The v e survival probability P ee from 
Eq. (|D reads then 

P ee = P & P (B + (l-P Q )(l-P e ) 



+2^P & (1 - P )P e (l - Pe) cose , (12) 
where the global oscillation phase is given by 

e=^(l-fe-*0-*e), (13) 



with the definitions 



3r ~ ~j~ ~ 4-7 x 10- s , (14) 

«• = ■ (16) 



4 Although Rq is relatively small (Rq/L = 4.7 x 10 3 ) , it is explicitly kept for later purposes. The 
Earth radius can instead be safely neglected (Rq/L = 4.3 x IO -5 ). 
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A relevant extension of the 2v formula ([12]) is obtained for 3z/ oscillations, required to 
accommodate solar and atmospheric neutrino data Assuming a third mass eigenstate 
2/3 with m 2 = |m 2 — m\ 2 \ 3> 5m 2 , the 3z/ survival probability can be written as 



where is the (z/ e , i^ 3 ) mixing angle, and P 2 ^ is given by Eq. (|T2"D, provided that the electron 
density A e is replaced everywhere by c^A^ (see [^0] and refs. therein). Such replacement 



implies that the 3v case is not a simple mapping of the 2u case, and requires sperine 
calculations for any given value of <p. 

We conclude this section by recovering some familiar expressions for P ee , as special cases 
of Eq. (|12"D. The JS limit (L mat /L osc — > 0, with complete suppression of oscillations inside 
matter) corresponds to P ~ c 2 ~ P ffi and to negligible 5 Q , <5 e . Then, neglecting also <5r, 
one gets from ( |ì~2]) the standard "vacuum oscillation formula," 



3>JS 

ee 



1 - sin 2 2o;sin 2 (7rL/L osc ) . (18) 



In the MSW limit (L/L osc — > oo), the global oscillation phase £ is very large and 
cos£ ~ on average. Furthermore, assuming for P & a well-known approximation in 
terms of the "crossing" probability P c between mass eigenstates in matter [in our nota- 
tion, P & ~ sin 2 u;^P c + cos 2 u;^(l — P c ), with oj^ calculated at the production point], one 
gets from ([12]) and for daytime (P e = c 2 ) the so-called Parke's formula |22| . 



P%Z ^ 1 + (ì ~Pc) cos 2u cos 2ul . (19) 



III. NUMERICAL TECHNIQUES 

In general, numerical calculations of the v transition amplitudes must take into account 
the detailed N e profile along the neutrino trajectory, both in the Sun and in the Earth. 

Concerning the Sun, we take N e from ("year 2000" standard solar model). Figure [I] 
shows such N e profile as a function of the normalized radius r/R®, together with its expo- 
nential approximation |l|] N e — N® exp(— r/ro), with A° = 245 mol/cm 3 and tq = P©/10.54. 
For the exponential density profile, the neutrino evolution equations can be solved analyt- 
ically [p|-[TT[|. In order to calculate the relevant probability P and the phase £ Q , we have 
developed several computer programs which evolve numerically the familiar MSW neutrino 
evolution equations along the Sun radius, for generic production points, and for any given 
value of Sm 2 /E £ [10~ 10 , 10~ 7 ] eV 2 /MeV and of tan 2 u. We estimate a numerical (fractional) 
accuracy of our results better than IO -4 , as derived by several independent checks. As a 
first test, we integrate numerically the MSW equations both in their usuai complex form 
(2 real + 2 imaginary components) and in their Bloch form involving three real amplitudes 
|24| , obtaining the same results. We have then repeated the calculations with different inte- 
gration routines taken from several computer libraries, and found no significant differences 
among the outputs. We have optionally considered, besides the exact N e profile, also the 
exponential profile, which allows a further comparison of the numerical integration of the 
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MSW equations with their analytical solutions, as worked out in |10|JTT[] in terms of hyper- 
geometric functions (that we have implemented in an independent code). Also in this case, 
no difference is found between the output of the different codes. 

Concerning the calculation of the quantities and in the Earth, we evolve analyti- 
cally the MSW equations at any given nadir angle rj, using the technique described in |25f , 



which is based on a five-step biquadratic approximation of the density profile from the Pre- 
liminary Reference Earth Model (PREM) |25j and on a first-order perturbative expansion of 



the neutrino evolution operator. Such analytical technique provides results very close to a 
full numerical evolution of the neutrino amplitudes, the differences being smaller than those 
induced by uncertainties in N e f25 |. In particular, we have checked that, for Sm 2 /E < IO -7 
eV 2 /MeV, such differences are ^5 IO -3 . In conclusion, we are confident in the accuracy of 
our results, which are discussed in the following sections. 

IV. MATTER EFFECTS IN THE SUN 

Figure ^ shows, in the mass-mixing piane and for standard solar model density, isolines 
of the difference c 2 — P Q (solid curves), which becomes zero in the just-so oscillation limit 
of very small 5m 2 /E. The isolines shape reminds the "lower corner" of the more familiar 
MSW triangle [P^| . Also shown are isolines of Constant resonance radius R Tes / R Q (dotted 
curves), defined by the MSW resonance condition L osc /L mat (-R ros ) = cos2u;. The values of 
c 2 — P & are already sizable (a few percent) at 5m 2 / E ~ 10~ 9 , and increase for increasing 
Sm 2 /E and for large mixing [tan 2 aj ~ 0(1)], especially in the first octant, where the MSW 
resonance can occur. The difference between matter effects in the first and in the second 



octant can lead to observable modifications of the allowed regions in fits to the data |L9 
and to a possible discrimination between the cases uj < | and uj > ? fll7| , [i9 



In the whole parameter range of Fig. |2], it turns out that, within the region of v production 
(t/Rq ^5 0.3), it is L ma t(r) <C L osc (and thus sin 2 2c<j^ ~ O).0 As a consequence, ali the 
curves of Fig. do not depend on the specific v production point (as we have also checked 
numerically) , and no smearing over the v source distribution is needed in the quasi- vacuum 
regime. This is a considerable simplification with respect to the MSW regime, which involves 
higher values of Sm 2 /E and thus shorter (resonance) radii, which are sensitive to the detailed 
v source distribution. 

Figure ^| shows, in the same coordinates of Fig. 0, the isolines of c 2 — P G corresponding to 
the exponential density profile (dotted curves), for which we have used the fully analytical 
results of p!0|JTT| . (Identical results are obtained by numerical integration.) The solid lines 
in Fig. refer to a well-known approximation (sometimes called "semianalytical" ) to such 
results, which is obtained in the limit N e — > at the Sun surface (it is not exactly so 
for the exponential profile, see Fig. [l]). More precisely, the zeroth order expansion of the 
hypergeometric functions fliQ| , |TTf in terms of the small parameter z = i\^2G FroN e (R Q ) (\z\ 



0.16) gives, for N e (0) — » oo, the result P & ~ [exp(7c 2 ) — l]/[exp(7) — 1] with 7 = TTr 5m 2 /E 
|IIJ , which leads to the QV prescription discussed in |19[ . The differences between the solid 



5 This is also indicated by the fact that R res /Ro ^ 0.55 in the Sm 2 /E range of Fig. 
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and dotted curves in Fig. |3| are essentially due to the "solar border approximation" (N e — > 0) 
assumed in the semianalytical case; indeed, the differences would practically vanish if the 
exponential density profile, and thus the "effective" Sun radius, were unphysically continued 
for r ^> R® (not shown). Such limitations of the semianalytical approximation have been 
qualitatively suggested by the authors of |27j but, contrary to their claim, our Fig. [3] shows 



explicitly that the semianalytical calculation of P® represents a reasonable approximation 
to the analytical one for uj in both octants, as also verified in 



A comparison of the results of Fig. ^| (true density) and of Fig. |] (exponential density) 
shows that, in the latter case, the correction term c 2 — P® tends to be somewhat over- 
estimated, in particular when the semianalytical approximation is used. We have verified 
that such bias is dominantly due to the difference (up to a factor of ~ 2) between the 
true density profile and its exponential approximation around r/R® ~ 0.8 (see Fig. [I]) and, 
sub dominantly, to the details of the density profile shape at the border (r/R® — » 1). As 
a consequence, the "exponential profile" calculation of P® (either semianalytic |~7|,|I11 or 
analytic) tends to shift systematically the onset of solar matter effects to lower values of 
Sm 2 /E. For instance, at tan 2 u; ~ 1 (maximal mixing), the value c 2 — P® = 0.05 is reached 
at 8m 2 /E ~ 8 x IO -10 eV 2 /MeV for the true density, and at 5m 2 /E a factor of ~ 2 lower for 
the exponential density. In order to avoid artificially larger effects at low Sm 2 /E in neutrino 
data analyses, one should numerically calculate P Q with the true electron density profile. 
The difference between the numerical calculation and the semianalytic approximation is also 
briefly discussed in |L9| for 5m 2 /E < 10~ 8 eV 2 /MeV (where P c ~ P®). 



Concerning the phase factor 8®, we confimi earlier indications [^III] about its small- 
ness, in both cases of true and exponential density. In the latter case, the semianalytic 
approximation gives [10,12|, for 8m 2 / E — > 0, the a;-independent result 



8® ~ L- 1 {r [\n{V2G F N%) + j E ] - R Q } , (20) 

where 7_b(— 0.577) is the Euler Constant. Numerically, 8® ~ —5.39 x IO -4 , much smaller 
than 5r. As far as 6m 2 /E < IO -8 (IO -7 ) eV 2 /MeV, we find that 5 e differs from such 
value by less than 20% (50%), the exact difference depending on the value of tan 2 u (not 
shown). Therefore, S & is an order of magnitude smaller than 5r in the whole mass-mixing 
range considered. A similar behavior is found by using the true density profile {\5 & \ ^ IO -3 
everywhere). As for P & , we find also for 5 Q no significant dependence on the neutrino 
production point in the Sun core. 

In conclusion, we find only modest differences between the analytical calculations of 
P Q , based on the exponential density profile, and its semianalytical approximation. The 
difference with respect to the numerical calculation of P®, based on the standard solar 
model profile, is instead more pronounced, and leads to a factor of ~ 2 difference in the 
value of 8m 2 / E where QV effects start to be significant. Therefore, we recommend the use 
of the electron profile from the standard solar model, implying numerical calculations of 
P®-f\ Concerning the phase S®, it simply shifts the global oscillation phase £ in Eq. (|TJ) by 



6 The interested reader can obtain numerical tables of P®, calculated for the standard solar model 
density, upon request from the authors. 
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less than one perniili, and thus it can be safely neglected in ali current applications. 



V. MATTER EFFECTS IN THE EARTH 

Strong Earth matter effects typically emerge in the range where L osc ~ P ma t within 
the mantle (N e ~ 2 mol/cm 3 ) or the core (N e ~ 5 mol/cm 3 ), as well as in other ranges 
of mantle-core oscillation interference |29|| , globally corresponding to 5m 2 /E ~ 10 _7 -10 -6 
eV 2 /MeV. Therefore, only marginai effects are expected in the parameter range considered 
in this work, as confirmed by the results reported in Fig. |j. 

Figure |4] shows isolines of the quantity c 2 — P e , which becomes zero in the just-so 
oscillation limit of very small 5m 2 /E. The solid curves corresponds to a nadir angle rj = 0° 
(diametral crossing of neutrinos) and the dotted curves to 77 = 45° (crossing of mantle only). 
For other values of 77 (not shown), the quantity c 2 — P© has a comparable magnitude. In the 
current neutrino j argon, the Earth effect shown in Fig. 0] is operative in the lowermost part 



of the so-called "LOW" MSW solution [16fl to the solar neutrino problem, or, from another 
point of view, to the uppermost part of the vacuum solutions [fL5| . Concerning the phase 
correction 5® (not shown), it is found to be smaller than 1.5 x 10~ 5 in the whole mass-mixing 
piane considered, and thus can be safely neglected. 

In practical applications, the correction term c 2 — P ffi must be time-averaged. This poses, 
in principle, a tedious integration problem, since such correction appears, in Eq. (|Ì~2|) , both in 
the amplitude of the oscillating term (oc cos£) and in the remaining, non-oscillating term. 
While the integration over time can be transformed, for the non-oscillating term, into a 
more manageable integration over r] |[25||, this cannot be done for the oscillating term, which 



depends on time both through the prefactor yP e (l — P ffi ) and through the phase £ (via 
eccentricity effects). However, in region of the mass-mixing piane where the Earth effect 
is non-negligible (namely, where c 2 — P© ~ few %), it turns out that P ~ 0, so that the 
total amplitude of the oscillating term is small. Moreover, as it will be discussed in the next 
section, energy smearing effects strongly suppress the oscillation factor cos£ in the same 
region. Therefore, for Sm 2 /E > 10~ 8 eV 2 /MeV, the oscillating term in Eq. (|I^) is doubly 
suppressed, and one can safely consider only the non-oscillatory terms. 



VI. DAMPING OF THE INTERFERENCE TERM 

In the just-so regime, the interference factor cos£ in Eq. can lead to an observable 
modulation both in the energy and in the time distribution of solar neutrino events. This 
modulation gradually disappears as 5m 2 / E increases, as a consequence of several decoher- 
ence effects ((cos£) — > 0) which are typically dominated by energy smearing |$0|. The 

broader the neutrino spectrum, the lower the values of Sm 2 /E where decoherence becomes 
importanti therefore, it suffices to consider the narrowest spectra (the so-called Be neutrino 
"lines" [fn]| ) for illustration purposes. 



Let us consider a v "line" energy spectrum s(E), with JdEs(E) = 1 and (E) 
JdE s(E) E. The energy-averaged oscillating factor, 
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H™^)W™ cos (^)' (2i) 

can be written, in the narrow-width approximation (AE = E — (E) <C (E)), in terms of the 
Fourier transform of the spectrum, 

s(r) = JdEs{E)e iAET . (22) 

More precisely, 

^/-^>»(i;('-3) <*> 

= J "»(w (1 -' ) ) 1 (24) 

where £) = |s(r)|, 5 = t^ 1 (E}^ 1 args(r), and r = 5m 2 L/2(E} 2 . 

Using the s(i£) profile for the Be lines from PTf , we find that the phase correction 5 
is smaller than a few perniili in the region where the damping factor D is greater than a 
few percent. If 5 is neglected, the average oscillation factor can be simply written as the 
oscillation factor at the average energy, times a damping terni D: 

I (5m 2 L\\ ^ (5m 2 L\ 



A similar approach to smearing effects for neutrino lines was developed in ||13|| . 

Figure |5] shows the damping factor D for the two Be neutrino lines at (E) = 0.863 
and 0.386 MeV. The factor is negligible for 5m 2 /E > IO -8 eV 2 /MeV, implying that the 
oscillation pattern is completely smeared out in such range.[] The onset of smearing effects 
is shifted to even lower values of Sm 2 /E for continuous energy spectra. In conclusion, the 
integration over neutrino energy makes the average oscillation factor (cos£) always negligible 
in the range where Earth matter effects are important, even for the narrow Be neutrino lines. 
Such accidental simplification should not make one forget that, in general, the transition 
from coherent to incoherent oscillations is a complex phenomenon that deserves further 



studies (see, e.g., [^2| for a recent ab initio approach). 



Finally, we recali that integration over time acts as a further damping factor [jTT . At first 
order in the eccentricity e, one has that L(t) ~ L[l — ecosa(t)] and à(t) ~ l + 2ecosa, where 
a is the orbitai phase, and t G [0, 2tt} is the normalized time variable (t = at perihelion). 
Then the yearly average of cos£ [times the square-law factor L 2 /L 2 (t)} can be performed 
analytically f33|, giving 



1 /** L 2 (5m 2 L(t)\ T / 5m 2 L\ (5m 2 L\ 



7 This range corresponds approximately to L osc /L < w/(E), where w ~ 0(1 keV) is the line 
width. 
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where Jq is the Bessel function, acting as a further damping term for large values of its 
argument. 

Notice that the maximum fractional variation of the orbitai radius, (L max — L min )/L = 
2e = 3.34 x 10~ 2 , is ari order of magnitude larger than Ór = Rq/L = 4.7 x 10~ 3 which, in 
turn, is larger than the phase corrections 5© and 5©. Therefore, one can safely neglect 5r, 
5q and 5© in practical applications involving yearly (or even seasonal) averages, as we do 
in this work. However, for averages over shorter time intervals, such approximation might 
break down. In particular, Sr (<5 q ) might be comparable to the monthly (weekly) variations 
of the solar neutrino signal. The observability of such short-time variations is beyond the 
present sensitivity of real-time solar v experiments and would require, among other things, 
very high statistics and an extremely stable level of both the signal detection efficiency and 
of the background. If such difficult experimental goals will be reached in the future, some of 
the approximations discussed so far (and recollected in the next section) should be revisited 
and possibly improved. 



VII. PRACTICAL RECIPES 

We have seen in the previous sections that, as Sm 2 /E increases, the deviations of P© 
(and subsequently of P©) from the vacuum value c 2 become increasingly important. We 
have also seen that the phase corrections 5 & and 5© are smaller than 5r = Rq/L, which 
can in turn be neglected in present applications, so that one can practically take the usuai 
vacuum value for the oscillation phase, £ ~ 5m 2 L/2E. We think it useful to organize known 
and less known results through the following approximate expressions for the calculation of 
P ee , which are accurate to better than 3% with respect to the exact, general formula (p~2[ ) 
valid at any 5m 2 /E. 

For 5m 2 /E <5x 1CT 10 eV 2 /MeV, one can take P© ~ P© ~ c 2 , and obtain the just-so 
oscillation formula 

Pe!^ct + st + 2 S 2 c 2 cos^ (27) 

with £ = 5m 2 L/2E. For 5 x 1(T 10 < 5m 2 / E < 1(T 8 eV 2 /MeV, one can stili take P© ~ c 2 , 
but since P© ^ c 2 (quasi-vacuum regime) one has that 



- <£Po + 4(1 - Po) + 2s uCu) JP e (l - P©) cose , (28) 



where £ = óm 2 L/2E, and P© has to be calculated numerically (see also [10— 12| for earlier 
versions of the above equation). 

Finally, for 5m 2 /E > IO -8 eV 2 /MeV, also Earth matter effects are important (P© ^ c 2 ); 
however, this complication is balanced by the disappearance of the oscillating term (cos £ ~ 
0) due to unavoidable smearing effects, so that the usuai MSW regime is recovered: 



>MSW 



P©P© + (1-P©)(1-P©) . (29) 



Time averages are then relatively simple to implement. In Eqs. (p7[) and (2^), the yearly 



averages of cos£ can be performed analytically [Eq. (JZBP]. In Eq. (^), the nighttime average 
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of P® can be transformed into a more manageable integration over the nadir angle, both for 



yearly [25] and for seasonal [34] averages. 

To summarize, the above sequence of equations describes the passage from the regime 
of just-so to that of MSW oscillations, via quasi-vacuum oscillations. In the JS regime, 
oscillations are basically coherent and do not depend on the electron density in the Earth or 
in the Sun (N e — > oo). In the MSW regime, oscillations are basically incoherent (L — > oo) 
and, in general, depend on the detailed electron density profile of both the Sun and the 
Earth. In particular, in the MSW regime one has to take into account the interplay between 
the density profile and the neutrino source distribution profile. The intermediate QV regime 
is instead characterized by partially coherent oscillations (with increasing decoherence as 
8m 2 / E increases), and by a sensitivity to the electron density of the Sun (but not of the 
Earth). Such sensitivity is not as strong as in the MSW regime and, in particular, QV effects 
are independent from the specific v production point, which can be effectively taken at the 
Sun center. 

For the sake of completeness, we mention that, for high values of 5m? / E (^> 10~ 4 
eV 2 /MeV), corresponding to L osc <C L mat in the Sun, the sensitivity to matter effects is 
eventually lost both in the Sun and in the Earth (P ~ P ffi ~ c 2 ), and one reaches a fourth 
regime sometimes called of energy-averaged (EA) oscillations, which is totally incoherent 
and iV e -independent: 

pI k * 4 + 4 . (30) 

Such regime, which predicts an energy-independent suppression of the solar neutrino flux, 
seems to be disfavored (but perhaps not yet ruled out) by current experimental data on 
total neutrino rates. In conclusion, for óm 2 /E going from extremely low values to infinity, 
one can identify four rather different oscillation regimes, 

JS -> QV -> MSW -> EA , (31) 

each being characterized by specific properties and applicable approximations. Experiments 
stili have to teli us unambiguously which of them truly applies to solar neutrinos. 

Vili. THREE-FLAVOR OSCILLATION ANALYSIS 



As discussed in |19| , in the QV regime the 1v survival probability fl28|) is non-symmetric 
with respect to the operation u — > | — u, which instead holds for JS oscillations [Eq. (|2"7j)].p| 
Here we extend such observation to 2>v oscillations, under the hypothesis 5m? <C m 2 which 
leads to Eq. (|ì~7j) . The 2u case is recovered for = 0. 

Equation ( [P7D for P^ v preserves the (a)symmetry properties of P<i v under the replacement 
u — > |— u. Therefore, while P 3 J ^ [obtained from Eqs. (|P7|) and (p7|) 1 is symmetric with respect 
to the uj — | value, the expression of P^) 1 [obtained from Eqs. (|T7D and (P5|)1 is not. Such 



8 In the MSW regime, the mirror asymmetry of the first two octants was explicitly shown in [20 1. 
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properties become evident in the triangular representation of the solar 3is mixing parameter 



space discusseci in p0|,|33|, to which the reader is referred for further details. 

Figure H shows, in the triangular plot, isolines of P^ v (dotted lines) for 5m 2 / E dose to 
1.65 x IO -9 eV 2 /MeV, corresponding to about 100 oscillation cycles. More precisely, the 
six panels correspond to £ = 100 x 2n + A£, with A£ from to tt in steps of ir/ 5. The 
dotted isolines are asymmetric with respect to u = f , as is the QV correction P & — c 2 . For 
increasing s\ (upper part of the triangle), the symmetry tends to be restored, mainly because 
the asymmetric term c^,P 2j , in Eq. (|T7|) is suppressed by the prefactor ci and, marginally, 
because the effective electron density c 2 iV e in P 2u is also suppressed.^] Exact mirror symmetry 
at any <fi is restored only if the QV corrections are switched off (solid lines); in such case, 
Pj, v is a quadratic form in the coordinates and s 2 , and the isolines are conical curves. In 
any case, P ee becomes typically too low (too high) in the centrai region (in the corners) of 
the triangle where, as a consequence, one should expect only marginai solutions to the solar 
neutrino deficit. 

The asymmetry with respect to u = f also shows up in solar neutrino data fits |L9| . 



Figure ^reports the results of our global (rates + spectrum + day/night) three-flavor analysis 
in the mass-mixing range 5m 2 G [10~ n , 10~ 8 ] eV 2 and tan 2 u; G [IO -2 , IO 2 ], for several 
representative values of tan 2 0. We only show 99% C.L. contoursQ (V DF = 3) for the sake 
of clarity. The theoretical [35] and experimental p6H4T[] inputs, as well as the x 2 statistica! 
analysis fl42" |, are the same as in Ref. [|l6j (where MSW solutions were studied). Here, 
however, the range of Sm 2 is lower, in order to show the smooth transition from the MSW 
solutions to the QV and finally the JS ones as Sm 2 is decreased. In particular, the solutions 
shown in Fig. |7| represent the continuation, at low 5m 2 , of the LOW MSW solutions shown 
in Fig. 10 of |jrj| (panel by panel). As anticipated in the comments to Fig. the mirror 
asymmetry around tan 2 cu = 1 decreases for decreasing Sm 2 (JS regime); a little asymmetry 
is stili present even at 5m 2 ~ 10~ 10 eV 2 , where the gallium rates (sensitive to E as low 
as ~ 0.2 MeV) start to feel QV effects. In the region where QV effects are important, 
the solutions are typically shifted in the second octant (u <^ 7r/4), since the gallium rate is 
suppressed too much in the first octant (see also fl9fl). A similar drift was found for the LOW 



9 When the effective electron density c^N e is small (s 2 — ► 1), the statement that the quasi- vacuum 
probability P 2v does not depend on the production point in the core (see Sec. IV) is not strictly valid, 
and one should in principle consider also the v source profile. However, since P 2v is correspondingly 
suppressed by c^, the source smearing effect is numerically small, and one can effectively discard 
it (e.g., by taking the production point at r = at any s^). 

10 At 95% C.L. the solutions would be mostly located at 5m 2 < 1CT 9 eV 2 . 

11 A technical remark is in order. The minimum value of x 2 in the piane of Fig. [?] (Xmin = 31.8) is 
reached at (5m 2 /eV 2 , tan 2 lo, tan 2 0) = (4.4 x IO" 10 , 2.4, 0.1). For 5m 2 > IO -8 eV 2 (MSW regime), 
the minimum value is Xmin = ^7.0 [|l6| for the same input data. In order to match the results of 
Fig. in this work and of Fig. 10 in juj], we nave adopted the absolute minimum {x^ìn = 27.0) to 
draw the Ax 2 = 11.34 contours in Fig. ^. 
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MSW solution |T(|. At any Sm 2 , the asymmetry decreases at large values of (tan 2 > 1.5) 



which, however, are excluded by the combination of accelerator and reactor data |2T] , unless 
the second mass square difference m 2 turns out to be in the lower part of the sensitivity range 
of the CHOOZ experiment (m 2 ~ IO -3 eV 2 ). For = 0, the standard two-flavor case 



is recovered, and the results are comparable to those found in ||15|| . The Super-Kamiokande 
spectrum plays only a marginai role in generating the mirror asymmetry of the solutions 
in Fig. ^, since the modulation of QV effects in the energy domain is much weaker than 
the one generated by the oscillation phase £. We find that, at any given óm 2 < IO -8 eV 2 , 
the x 2 difference at symmetric lo values is less than ~ 1 for the 18-bin spectrum data fit. 
Therefore, QV effects are mainly probed by total neutrino rates at present. 

We think it is not particularly useful to discuss more detailed features of the current 
QV solutions, such as combinations of spectral data or rates only, fits with variations of hep 
neutrino ffux, etc. (which were instead given in |Tj| for MSW solutions). In fact, while the 
shapes of current MSW solutions are rather well-defmed, those of JS or QV solutions are 
stili very sensitive to small changes in the theoretical or experimental input. Therefore, a 
detailed analysis of the "fine structure" of the QV solutions in Fig. [7] seems unwarranted at 
present. 

Finally, in Fig. |8] we show sections of the allowed 3u solutions (at 99% C.L.) in the triangle 
representation, for six selected (increasing) values of Sm 2 . Solutions are absent or shrunk 
at s 2 ~ 0.5, where the theoretical v flux underestimates the gallium and water-Cherenkov 
data. The lowest value of 5m 2 (0.66 x 10~ 10 eV 2 ) falls in the JS regime, so that the ring-like 
allowed region (which resembles the curves of iso-P ee in Fig. is symmetric with respect to 
the vertical axis at lo = 7r/4. However, as 5m 2 increases and QV effects become operative, 
the solutions become more and more asymmetric, and shifted towards the second octant of 

LO |7|,0. 



Figures [7| and § in this work, as well as Fig. 10 in [Tj| , show that solar neutrino data, by 



themselves, put only a weak upper bound on the mixing angle 0. Much tighter constraints 
are set by reactor data []43"| |, unless the second mass square difference m 2 happens to be 
< 10~ 3 eV 2 (which seems an unlikely possibility). In any case, QV effects are operative also 
for a small (or zero) value of sì. 



IX. CONCLUSIONS 

We have presented a thorough analysis of solar neutrino oscillations in the "quasi- 
vacuum" oscillation regime, intermediate between the familiar just-so and MSW regimes. 
The QV regime is increasingly affected by matter effects for increasing values of Sm 2 . We 
have calculated such effects both in the Sun and in the Earth, and discussed the accuracy 
of various possible approximations. We have implemented the QV oscillation probability in 
a full three-flavor analysis of solar neutrino data, obtaining solutions which smoothly join 
(at 5m 2 ~ IO -8 eV 2 ) the LOW MSW regions found in Jl6| for the same input data. The 



asymmetry of QV effects makes such solutions different for lo < | and lo > |, the two cases 
being symmetrized only in the just-so oscillation limit of small 5m 2 . 
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FIGURES 



FIG. 1. Radiai profìle of the electron density in the Sun from the standard solar model (solid 
line), together with its exponential approximation (dashed line). 

FIG. 2. Isolines of the solar correction terni c 2 — Pq in the mass-mixing piane (solid curves), 
for standard solar model density ^3|. Isolines of MSW resonance radii (dashed curves) are also 
shown. 

FIG. 3. As in Fig. |2|, but for the exponential density profile. Dotted and solid lines correspond 
to analytical calculations and to their semianalytical approximation, respectively. 

FIG. 4. Isolines of the Earth correction term c 2 — P® in the mass-mixing piane, for PREM |26| 
density. The solid and dotted lines correspond to a nadir angle equal to 0° and 45°, respectively. 

FIG. 5. Damping factor for the oscillating term, due to the finite Be line width. See the text 
for details. 

FIG. 6. Three-flavor oscillations: Isolines of P ee in the triangular representation of the solar v 
mixing parameter space, for representative values of the oscillation phase £. The dotted and solid 
lines refer, respectively, to calculations with and without quasi-vacuum effects. When such effects 
are included, the mirror symmetry tu — > ^ — lo is broken. 

FIG. 7. Three-flavor oscillations: Contours at 99% C.L., as derived from a global analysis 
of neutrino data for 5m? < IO -8 eV 2 (quasi-vacuum regime). The asymmetry of the solutions is 
evident for increasing values of 5m 2 . The solutions represent the continuation (at low Sm?) of the 
MSW 99% C.L. regions reported in Fig. 10 of |§. 

FIG. 8. Sections of the 3v allowed volume (99% C.L.) for six representative values of Sm 2 
(0.66, 2.36, 6.42, 8.76, 10.8, and 29.5 x IO" 10 eV 2 ), shown in the solar triangle plot. 
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FIG. I]. Radiai profile of the electron density in the Sun from the standard solar model (solid 
line), together with its exponential approximation (dashed line). 
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FIG. ^. Isolines of the solar correction terni — Pq in the mass-mixing piane (solid curves), 
for standard solar model density p3|. Isolines of MSW resonance radii (dashed curves) are also 
shown. 
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FIG. |||. As in Fig. ^, but for the exponential density profile. Dotted and solid lines correspond 
to analytical calculations and to their semianalytical approximation, respectively. 
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ci — Rg for PREM electron density 
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FIG. ^. Isolines of the Earth correction term — in the mass-mixing piane, for PREM |26| 
density. The solid and dotted lines correspond to a nadir angle equal to 0° and 45°, respectively. 
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Damping factors for Be lines 
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FIG. ||. Damping factor for the oscillating terni, due to the finite Be line width. See the text 
for details. 
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3v mixing: curves of iso-P(z/ e — > v s ) 
for oscillation phase = 10O27T + A£ 




FIG. ||. Three-flavor oscillations: Isolines of P ee in the triangular representation of the solar v 
mixing parameter space, for representative values of the oscillation phase £. The dotted and solid 
lines refer, respectively, to calculations with and without quasi-vacuum effects. When such effects 
are included, the mirror symmetry uj — > f — lo is broken. 
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FIG. 0. Three-flavor oscillations: Contours at 99% C.L., as derived from a global analysis of 
neutrino data for 5m 2 < IO" 8 eV 2 (quasi- vacuum regime). The asymmetry of the solutions is 
evident for increasing values of óm 2 . The solutions represent the continuation (at low dm?) of the 
MSW 99% C.L. regions reported in Fig. 10 of p|. 
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Solar 3v analysis (99 % C.L.) 
óm 2 /10~ 10 eV 2 = 




FIG. ||. Sections of the 2>v allowed volume (99% C.L.) for six representative values of dm 2 (0.66, 
2.36, 6.42, 8.76, 10.8, and 29.5 x 10~ 10 eV 2 ), shown in the solar triangle plot. 
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